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SUMMARY 


A detailed exploration of the field of mean and fluctuating quan- 
tities in a two-dimensional turbulent channel flow is presented. The 
measurements were repeated at threie Reynolds numbers, 12,300, 30,800, 
and 61 , 600 , based on the half width of the channel end the may-lmum 
velocity. A channel of 5-inch width and 12:1 aspect ratio was used for 
the investigation. 

Mean— speed and axial— fluctuation measurements were made well within 
the laminar sublayer. The . semi theoretical predictions concerning the 
extent of the laminar sublayer were confirmed. The distribution of the 
velocity fluctuations in the direction of mean flow u' shows that the 
influence of the viscosity extends farther from the wall than indicated 
by the mean velocity profile, the region of influence being approxi- 
mately four times as wide. 


.Fluctuations perpendicular to the flow in the lateral and vertical 


directions 


cient k = 


u' v 


and w' , respectively, and the correlation coeffi— 

were also measured. The turbulent shearing stress 


was computed by three different methods.' The results show satisfactory 
agreement for the two lower Reynolds numbers. In the case of the high est 
Reynolds number, however, the total shearing stress r obtained from the 
fluctuation measurements was approximately 20 percent lower than tha t 
computed from the mean— velocity and mean— pressure measurements. All 
dimensionless mean fluctuating quantities were found to decrease with 
increasing Reynolds number. Measurements of the scales of turbulence Iy 
and Lj and microscales' of turbulence Xy and X z across the channel 
are presented and their variation with Reynolds number is discussed. 

Using a new technique, values for the microscale X^ were obt ain ed; a 
new method for estimating the scale 1^. is also given. 
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The energy balance in the turbulent flow field was calculated from 
the measured quantities. From this calculation it is possible to give 
a descriptive picture of turbulent-energy diffusion in the center 
portion of the channel cross section. 

For the flow corresponding to a Reynolds number of 30,800 the 

energy spectrum of the u' 2 — fluctuations at various points across the 
channel, including one in the laminar sublayer, was obtained. 

At station ^ = 0.4, where y is lateral distance and d is half 
d 

width of the channel, the contribution to the turbulent shear stress 
from various frequency bands was measured and it was found that the con- 
tribution corresponding to frequenci es above 1500 cycles per second is 

negligible. Since the spectrum of u‘ 2 at this point extends to about 
5000 cycles per second it is evident that the high frequencies are 
nearly isotropic in agreement with Kolmogorof f * s hypothesis. 


INTRODUCTION 


In recent years a considerable step forward was made in the theory 
of turbulence. Kolmogoroff (reference l), Heisenberg (reference 2), 
and Onsager (reference 3) obtained independently an energy-spectrum law 
that holds in rather restricted types of turbulent flows. This progress 
gave a new impetus to both -theoretical and experimental investigations. 
The experimental worker may follow two principal methods of approach to 
the problem. First, he may establish flow fields which satisfy suffi- 
ciently the assumptions of the new theory, namely, that the flow is 
isotropic and of high Reynolds number so that the influence of viscosity 
■ is a minimum and the effect of the turbulence-producing mechanism is 
small. Under these conditions he may measure quantities such as corre- 
lation functions, scales, and microscales that axe defined exactly in 
the flow field and may compare his results with those predicted by the 
theory. The main difficulty with this method is to predict how closely 
one has to appro xima te in the 'actual flow the conditions assumed in the 
theory. In other words, the sensitivity of the theory to deviations 
from true isotropic conditions is not known and in case the measurements 
do not agree with the theory one does not know whether to attribute the 
disagreement to faulty assumptions in the theory or to the incomplete 
isotropic conditions of the flow. Furthermore, the imposition of the 
condition of isotropy to fluctuating fields may be a very strong 
restriction and the study of such fields may not yield the complete 
picture of the turbulence mechanism. 
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The second method of approach is to establish a simple nonisotropic 
turbulent flow field with well-defined boundary conditions and, in the 
light of the existing theories, to try to obtain information on the 
mechanism of energy transfer from the low frequencies of the energy 
spectrum to the high ones. The present investigation of a turbulent 
channel flow has this purpose in mind in its long-range program. 

The difficulties of this method are Immediately realized. Because 
of the nonisotropic nature of the flow, characteristic quantities such 
as scale and microscale are no longer well defined enri it is very 
possible that in the first stages of the investigation certain quan- 
tities will be measured that later will prove to be trivial. On the 
other hand, the main advantage of a fully developed c hann el flow is 
the fact that, in contrast with the flow behind grids, flow conditions 
are steady; no decay of mean or fluctuating quantities exists in the 
direction of the flow. Consequently, the turbulent energy goes through 
all of its stages of transformation across the channel section — 
turbulent-energy production from the mean flow, energy diffusion, and 
turbulent and laminar dissipation - and one may study these transforma- 
tions in detail. 

It has become clear that the phenomenological theories of turbulence, 
such as the mixing-length theories., have lost most of their importance. ' 
These theories, developed in the late twenties and early thirties, wel-e 
aimed specifically at an evaluation of the mean-velocity distribution in 
turbulent flow. The existing experimental evidence shows clearly that ' 
the mean-velocity distribution is very insensitive to the essential 
assumptions introduced into the phenomenological theories. In fact, 
purely dimensional arguments generally suffic'e to give the shape of the 
mean velocity profile with sufficient accuracy. For the further develop- 
ment of an understanding of turbulence, detailed measurements of the 
field of fluctuating rather than of mean velocities are necessary. The 
program of experimental research of which this work is one part is based 
on this reasoning. It is quite apparent and natural that the same con- 
clusions have been drawn by workers in the turbulent field elsewhere 
and that, in general, the main emphasis of current experimental inves- 
tigation is the exploration of the field of the fluctuating— velocity 
components. 

The present set of experiments deals with flow in a two-dimensional 
channel, that is, with pressure flow between two flat walls. Channel 
flow of this type is the simplest type of turbulent flow near solid 
boundaries which can be produced experimentally. The simpler Couette 
type flow requires that one wall move with constant velocity, a con- 
dition which is difficult to realize experimentally. It can be 
approx ima ted by the flow between concentric cylinders , but complications 
due to centrifugal forces arise here. The simple geometry of a 
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two-dimensional c hann el allows an integration of the Reynolds equations , 
and the turbulent shearing stress can then be related directly to the 
shearing stress on the surfaces which, in turn, can be determined from 
the mean-pressure gradient or the slope of mean velocity profile at 
the wall. 

The relation between the apparent stresses and the wall shearing 
stress can be usfed to advantage in two fashions. It is here possible 
to obtain the magnitude of the correlation coefficient responsible for 
the apparent shear by measuring only the intensities of the turbulent 
fluctuations. The turbulent shearing stress can also be measured 
directly by means of the hot-wire anemometer. A comparison with the 
shear distribution obtained from mean-pressure-gradient or mean— 
velocity-profile measurements serves then as a very useful check on 
the underlying assumptions, on the one hand, and specifically as a 
check on the reliability and accuracy of the direct measurements, on 
the other. 

Measurements of channel flow have previously been made by Doench 
(reference 4) , Nikuradse (reference 5), Wattendorf and Kuethe (refer- 
ences 6 and 7), and Reichardt (references 8 and 9). A few unpublished 
measurements have been made recently at the Polytechnic Institute of 
Brooklyn. 

Doench' s and Nikuradse 1 s measurements were concerned only with 
the mean— velocity distribution and are thus of' not too much interest for 
a comparison with the present set of measurements. Wattendorf measured 
the intensity of the fluctuating-velocity components and then deduced 
the correlation coefficient from the mean— pressure measurements. The 
technique for measuring the axial component of the velocity fluctuation 
was well developed at the time, but the cross component was only 
tentatively measured. 

The most complete set of measurements is due to the work of 
Reichardt, who measured velocity fluctuations in the direction of the 
flow and normal to the wall as well as the turbulent shear directly. 
Reichardt found very good agreement between the shearing stress deter- 
mined in these two ways; his paper comes closest to the present inves- 
tigation and his results will be used for comparison. The Reynolds 
number in Reichardt' s measurements was 8000,' which is lower than the 
range covered by the present measurements. A criticism which can be 
made of Reichardt' s investigation is his use of a tunnel of only 
1:4 aspect ratio. The two-dimensional character of the flow is thus 
somewhat doubtful. Wattendorf* s experiments were made in a channel of 
very large aspect ratio (l8:l) and are thus free of’ this criticism. 

In a preliminary investigation a channel 1 inch wide and 60 inches 
high was chosen. The measurements, however, have shown that in this 
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case the scale of turbulence is so small that great care must be taken 
to correct the hot-wire readings for the effect of wire length. Measure- 
ments of the microscales were, in fact, impossible in the 1— inch channel. 
The microscales were about 0.1 centimeter and thus smaller than the 
length of the wire. The corrections for the case of measurements of v* 
and w' were about 30 percent. Since the method of correction becomes 
very Ina ccurate for such large ratios of wire length to microscale, the 
measurements were repeated in a 5— inch channel with a 12:1 aspect ratio. 
This ratio was still large enough to insure two-dimensional flow' and the 
length corrections were greatly reduced. 

This investigation was conducted at the Guggenheim Aeronautical 
Laboratory, California Institute of Technology, under the sponsorship 
and with the financial assistance of the National Advisory Committee 
for Aeronautics. 

The author wishes to acknowledge the constant advice and help 
in both experimental methods and interpretation of results of 
Dr. H. W. Liepmann during this investigation. He also wishes to thank 
Dr. C. B. Millikan for his continuous interest in this research. The 
cooperation' of Dr. F. E. Marble and Mr. F. K. Chuang is much appreciated. 


SYMBOLS 


x longitudinal coordinate in direction of flow; x = 0 corre- 

sponds to channel exit 

y lateral coordinate; y = 0 corresponds to channel wall 

z vertical coordinate 

d half width of channel (2.5 in. - ) 

5 thickness of laminar sublayer 

u mean velocity at any point in channel 

U 0 maximum value of mean velocity 

u 1 instantaneous value of velocity fluctuations in direction of 

mean flow x 

v* ,w* instantaneous values of velocity fluctuations normal to mean 

flow in directions y and z, respectively 
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u' 


V 1 yW* 


P 

P 

T 


U T 

d 

V 

E 


root-mean— square value of velocity fluctuations in 
direction of mean flow x 

root-mean-square values of velocity fluctuations normal 
to mean flow in directions y and z, respectively 

pressure at any point in channel 
air density 
totad shearing stress 
shearing stress at wall 


ii; y (\/ T o/p) 


friction velocity 

absolute viscosity of air 
kinematic viscosity 

Eeynolds number based on half width of channel and maximum 
mean velocity 


correlation coefficient responsible for apparent 


shear 



E^EyjEjj. correlation coefficients as functions of X, Y, and Z, 
respectively 

microscales of turbulence 


u |g _ / du'\ 
XJ2 " W / ’ 


u'£ _ i/a u .\2 u ,2 

2 ■ Z ^ y ~) ’ 1 2 


V' 


= lfau'V 

“fez 
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hu’^y’h. 


scales of turbulence 





X,Y,Z distances (in the x— , y— , and z-direction, respectively) 

between points at which correlation fluctuations are 
measured 


0 


ratio of square of compensated and uncompensated 
fluctuations 


M 

n 


time constant expressing thermal lag of hot-wire 
frequency , cycles per second 


F^t 2 (n)>F(n) fraction of turbulent energy u’ 2 associated with band 

width dn 


F — — j (n) fraction of turbulent shear u'v* associated with band 

U width dn 


W dissipation of turbulent energy 

t time 


e lj®2 


voltage fluctuations 


ANALYTICAL CONSIDERATIONS 

Equations of Motion for Two-Dimensional Channel Flow 


The mean— and fluctuating— velocity components are denoted by uf 
and ui', respectively; Cartesian coordinates are designated by x^; 
and Pik denotes the components of the stress tensor which includes 
both the viscous and apparent (Reynolds) stresses. The Reynolds equa- 
tion and the continuity equation for steady flow can thus be written, 
using Cartesian tensor notation, 


P% 


du^ 


dPik 

dxk 


(la) 





(lh) 



where 
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where x, y, and z replace xi, and u’, v’, and w' replace \i i t 

for convenience in writing the now much simpler equations. Equation (lb) 
is thus automatically satisfied and equation (la) becomes: 



In fully developed turbulent flow the variation of mean values of the 
fluctuating quantities with x should be zero; that is, the flow 
pattern is independent of the streamwise direction. Hence,, equations (2) 
become: 


or 


1 dp _ d^u _ du* v r 

P dx ” ay2 dy 


1 <3p _ 1 dr 

P dx P dy 


(3a) 


1 dp dv 1 ^ 

P dy dy 


(3b) 


Differentiating equation (3b) with respect to x gives 


1 d 2 p 

p dxdy 


consequently dp/dx is independent of y and equation (3a) is immedi- 
ately integrable. Thus 


y dp du -j— . . _ ^ J 

i-_£. = v u'v' + Constant 

P dx dy 
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In the center of the channel the shear vanishes; hence if the channel 
has the width 2d, 

Constant = — 

P .dx 


and hence 


d) = v — u'v* 


i^/y _ 
P d^ 7 


Or in nondimens ional form: 

y - d dp v <i u u'v' 

pU 0 2 dx u 0 2 ^ U 0 2 

In terms of the shearing stress at the wall. 


( 4 ) 


to y — d _ v du u 


1 T7-t 


PUo 


U„ 2 dy U 2 


(5) 


It is evident from equations (4) and (5) that t q can he determined in 
three different ways: 

(a) From the mean-pressure gradient 

dp 


t d = -d 


dx 


(b) From the slope of the mean velocity profile near the wall 


to 


/du\ 

W/ y= o 


(c) From a direct measurement of u'v' 


To 


= ^pu'v' 


du \ 

— d — 1 


dy/y - d 
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since if y/d is not too small. 


M- — « pu 1 V 1 


and item (c) "becomes 


pu* v 1 




The technically simplest way to determine T Qy and thus in fact the 
complete shear distribution, is a measurement of bp/bx. This is the 
method applied in most investigations. 

To determine the slope of the velocity profile near the vail, the 
profile has to be known up to points very close to the wall, that is, 

at least to distances — = 10“^. This. in general requires the use of 

d 

the hot-wire anemometer and very precise measurements* 

The third method requires a direct measurement of the correlation 
between the axial and lateral velocity fluctuations. The technique of 
this type of measurement is known and was first applied by Reichardt 
and by Skramstad and, in somewhat different form, in recent investi- 
gations at the National Bureau of Standards, Polytechnic Institute of 
Brooklyn, and California Institute of Technology. 

In the present Investigation all three methods. have been applied 
and the results compared, with the exception of the flow at the lowest 
Reynolds number. In this case the pressure gradient is extremely small 
(approx. 0.0003 111131 of alcohol/cm) and reasonably accurate measurements 
were not possible. This comparison of the three methods has the 
advantage that it gives a good indication of the absolute accuracy of 
measurements of the fluctuating quantities and the correlation coef- 
ficient k. 


Energy Equation 

Writing equation (la) in the form 

pu k + M . v 2 Ui 

OXj£ ox^ 


( 6 ) 



12 


NACA TN 2123 


end multiplying it by up* the following relation is obtained: 

- P — v - t — - + uui ^ 

2 ox k oxp oxjOXj 

Transforming the last term by partial integration the energy equation' 
of the mean flow becomes 

1 ^i u i u k = _ ^i + 1 ^ 2u i u i _ /^i\p u i\ 

2 P dx k cfccp + 2 P Sxj^xj ^^XjJ^xJ ' 

The velocity perturbations can now be introduced: 

up = u + u' 

I 

U2 = v’ 


u 3 = w ' 


Since the velocities are independent of the coordinates x and z the 
following equation is obtained after averaging: 


-j— r du du’v' 1 d ,/,2 .2 ,o\ dp dv*P 

dy dy 2 dy x • dx dy 


d 2 u 1 d 2 7~TZ io 7o\ 

uu — - + — |i (u 1 ^ + v’c + — |il^ IItt 

"" " ~ \ ' \OXA /\OXj 

\ ^ / \ J J 


dy 2 2 dy 2 
Making use of equation (3a) this simplifies to 


du d 
dy dy 


pv' 


^ U* 2 + V* 2 + w* 2 ^ 


+ v'p 




+ v’ 2 + w* 2 ) 



NACA TN 2123 


13 


This form was obtained by Von Karm&n (reference 10 ) while discussing a 
nonisotropic flow in terms of the statistical theory of turbulence. In 
order to see the relative orders of magnitudes of the different terms 
in the equation. Von K&rtn^n expresses them in terms of a single 

velocity q = \Ju' 2 + v' 2 + w' 2 and a characteristic length D corre- 
sponding to the width of the channel. Since 

^ ta -u'v' = 0(q2) 


and 



the above equation may be written as 


q3 q3 
Ai — = An — + 
D 2 D 


u + A, v£ + . Bi 

3 d 2 V ^ X 2 


where the coefficients have the characteristics of correlation functions 
and X is the microscale of turbulence. For a high Reynolds number 

flow — » 1, it follows that Thus the second and third 

v D2 D 

term on the right side of the equation may be neglected. Since the 

pertinent quantities have been measured lurin'- the present work, the 

order of magnitude of these terms can be directly evaluated ami the 

omission of the terms is found to be justified. The above equation thus 

contains only terms of the forma q3/D and vq 2 A 2 ; these terms should 

be of the same order of magnitude and therefore = 0(l). It is 

vD 

of considerable interest to see whether experimental results confirm 
this relation. Taking as an example results from the measurements 

/ 3 2 

at R = 30,800 and 2. - 0.5 (where y « 11 x lO^t and ** 1.7 x H)3 

q = 52 centimeters per second 
Xy = 0.5 centimeter 


therefore 

qx£ 52 X 0.25 
VD 0.155 X 12.7 
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This ratio seems to be fairly constant for different Reynolds numbers 
and across the channel cross section. 

In view of the above dimensional considerations, the energy- 
equation may be written 



The equation expresses the fact that the energy produced by the 
turbulent shear forces at a certain point is partly diffused and partly 
dissipated. This equation has been used in estimating the energy 
diffusion across the channel, since from the measured quantities the 
production term and dissipation can be calculated. 


EQUIPMENT AND PROCEDURE 
Wind Tunnel 


The investigation was carried out in the wind tunnel shown in 
figure 1. The turbulence level is controlled by a honeycomb and seam- 1 
less precision screens, followed by an approximately 29:1 contraction. 
The screens have 18 meshes per inch and a wire diameter of 0.018 inch. 
The honeycomb consists of paper mailing tubes, 6 inches long and 1 inch 
in diameter . 


The over— all length of the channel is 23 feet. At the entrance 
section it is 3 inches wide and has an aspect ratio of 20:1. In a 
distance of 7 feet it expands to a width of 5 inches and the aspect 
ratio is reduced to 12:1. The walls of the exit portion of the channel 


o 

(about 6 ft) are made of jj- —inch— thick plywood with a ——inch birch 

inside cover. They "were specially treated in order to acquire a smooth 
finish and were reinforced in order to avoid warping (fig. 2); in spite 
of this a few percent of width variation existed. 


The tunnel is operated by a 62-horsepower stationary natural— gas 
engine — normally operating at a fraction of its rating — which drives 
two eight-bladed fans. The speed is remotely controlled by means of a 
small electric motor which drives the throttle through a gear and lead- 
screw system. 


The experiments were carried out at speeds of 3, 7*5, and 15 meters 
per second. 
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Traversing Mechanism 

Figure 3 shows the type of traversing mechanism used during the 
experiments. It consists simply of a micrometer screw on which the 
hot-wire support is fastened. The support can rotate freely in a plane 
perpendicular to the air flow; thus, the hot-wire may be adjusted 
exactly parallel to the wall. 

The zero reading of the traversing mechanism (y = 0) was care- 
fully found using the following method: The hot-wire was placed close 

to the wall (approx. 0.025 cm away); the distance between the wire 
and its image in the polished wall was measured by means of an ocular 
micrometer. The position of the wire is, of course, one-half of the 
observed distance and could be determined with an accuracy of. 

±0.0005 centimeter. 

In order to obtain the pressure distribution along the middle of 
the channel a 6-inch— long thin-walled tube 3/8 inch in diameter was used 
with a small static-pressure hole drilled close to its end. The tube 
was free to slide in a support at the entrance of the channel so that 
the position of the pressure hole relative to the channel exit could be 
changed. The tube was kept straight and under tension by means of 
weights. - 


Hot-Wire Equipment 

All velocity measurements were made with hot-wire anemometers. 

The frequency response of the amplifier— compensator unit of the hot- 
wire equipment, using a 0.0002^— inch wire at standard operating qon— 
ditions, is flat from approximately 2 to 10,000 cycles per second. 

The compensation of the hot-wire for thermal lag is accomplished 
by a capacity network. The range of time constants was chosen from 0 
to 1 millisecond corresponding to the characteristics of 0.00005— to 
0.00021*— inch wires at the operating conditions employed in general at 
GALCIT. No attempt was made to extend the range of compensation to 
larger values of M, since in this case the noise level soon becomes 
appreciable . * 

The correct' setting of the compensating unit was found using the 
square-wave method described by Kovas.znay (reference ll). The t ime 
constants of the wires used for turbulence measurements fell between 0.1 
and 0.9 millisecond. 

The output readings were taken with a thermocouple and precision 
potentiometer. 
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Mean-speed, measurements .- A "mil (0. 0005— in.) platinum wire of 

1— centimeter length was used, to measure the mean speed.. The measure- 
ments were made by the constant— resistance method. This method has the 
advantage of keeping the wire temperature constant through the velocity 
field. 

Turbulence measurements .— For the investigation of turbulent fluc- 
tuations, 0 . 00024— inch Wollaston wire was used. The wire was soft- 
soldered to the tips of fine sewing needles after the silver coating 
ha d been etched off. The two lateral components v' and w* and the 

correlation coefficient k = u T v* — were measured using the X-wire 

technique. The method was essentially the same as described in detail 
in previous reports from this laboratory. (See, for example, refer- 
ence 12.) 

The x-meters for shear measurements had angles between the wires 
of approximately 90°. The angles of the v',w' -meter were of the order 
of 30°. The wire length of the u* -meter was about 1,5 millimeters and 
that of the v' -meter was 3 millimeters; 

The parallel-wire technique used recently at the Polytechnic 
Institute of Brooklyn was also tried in order to obtain the lateral 
component of the velocity fluctuation and the shear close to the wall . 
The parallel-wire instrument is superior to the x-meter since it allows 
an exploration of the turbulent field up to distances of a few thous- 
andths of an inch from the wall. However, it was found impossible to 
obtain reliable values with this instrument. The corrections due to 
u'— fluctuations and to unequal heating of the wires were pronounced and 
not easily accountable. The method was therefore temporarily abandoned 
after considerable time and labor had been spent on it. 

Me asurements of correlations between two points .— The correlation 

u* at points along the y— and z-axis, 

u* (O)u' (Y) 

Ju’ (0) 2 ^u' (Y) 2 


u* (o)u' (z) 


functions between values of 
that is , 
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were measured using the standard technique (reference 12). The scales of 
turbulence Iy and L z and the microscales Xy and X, z at different 

points across the channel were obtained from these measurements. 

Measurements of X. x .— A method suggested by Townsend (reference 13) 
was applied to measure \ x . Using an electronic differentiation circuit 
the amplffier output signal was differentiated and X, x computed from 
the following relation: 



(■} r) 

The error involved in the approximation ^ > u 2_ ca n be esti— 

ot ox 

mated and it may be seen that the above relation holds with reasonable 
accuracy over the large center portion of the channel. 

A new technique for the measurement of X, x has also been applied. 
The method is described in detail in reference l4 and consists in 
counting the zeros of an oscillograph trace of the u'— fluctuations. 

From these counts \ x may be calculated directly by assuming a normal 
and independent distribution for both u' and du’/cfac: 

— = — X Average number of zeros of u* per second 
^x Uo 


It is known that the distribution of u' is closely a Gaussian one 
even in nonistropic turbulence (see, for instance, reference 15); 
however, for the case of du' fdx a small deviation from the normal dis- 
tribution was found (reference 13)* For the preliminary measurements 
of X, x reported presently, no corrections were applied as yet for this 
effect. Figure 4 shows an oscillograph trace of the u*— fluctuation in 
the middle of the channel at B = 30,800. The trace represents an 
interval of approximately l/20 second. 

Measurement of Lx.— The following simple procedure-*- was applied 

to obtain a rough estimate of scale of turbulence corresponding to 
correlations between points along the x-axis: Denote by F(n) the 


■^This method was suggested by Dr. H. W. Liepmann. 
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fraction of turbulent intensity which is contributed by 'frequencies 
between n and n + dn; that is. 


u* (n) 2 dn = u ,2 F(n) dn 


and thus 



dn = 1 


Consider now an uncompensated hot-^wire. If the time constant of this 
uncompensated wire is M, the response will be 



uncomp . 


where 



The total intensity for the uncompensated wire will then be given by 


2 



uncomp. 



F(n) dn 
1 + M 2 n 2 


Thus, if the ratio of the response of the same wire compensated and 
uncompensated is denoted by a. 


1 

CT 


r *<»> 

'-'0 1 + r ^ n 2 


(9) 


In order to estimate 


Lx, F(n) is assumed to have the simple form 


F(n) 


2 

7t 


u 0 


i 


+ n 


2 hL. 

Uo 2 


formulas of this general type have been proposed by Kamp^ de Feriet 
axiA by Frehkiel for determining the spectrum of turbulence from uncom- 
pensated hot-wire measurements by varying M (reference 1 6 ) . 
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which was given- by Dryden (reference 17) and corresponds to an experi- 
mental correlation curve. 


Then 


or with 


then 


Hence 



T) = 



a 


^x. 


1 

a 


p 00 


2 

rt 




Uo 


(l + T! 2 )(l + a 2 ri 2 ) 


1 

1 + a 


kjc - 


a - 1 


Thus by measuring the ratio of the mean squares of the fluctuating 
velocities with and without compensation, 1^ can be estimated. 

The procedure was checked in the flow behind a i-- inch grid where 

the lateral scale was measured and where because of isotropy the rela- 
tion = 2Iy should hold fairly well. 

The results were: 

a = 1.85 

M = 7 X 10 - ^ second 


Uo' = 1.5 X 10 3 centimeters per second 
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therefore 


I>x = 1.24 centimeters 

Direct measurements gave Ly = 0.553 centimeter. Hence 
Lx = 1.11 centimeters which shows that this method of estimating the 
value of Lx is satisfactory. 

Measurement of turbulent-energy spectrum .— In order to be able to 
measure accurately the energy distribution over a wide frequency band, 
it was found necessary to use extremely thin wires with small time 
constants and thus to increase the signal— to— noise 'ratio . Wires of 
0.00005— inch diameter, operating with time constants of approxi- 
mately 2 X 10“^ second, were used. In this way it was possible to 
detect energy values in the high-frequency bands that were 10~ 7 times 
the value corresponding to zero frequency. 

The hot-wire signal was fed into a. Hewlett-Packard wave analyzer, 
the output circuit of which was somewhat altered in order to be able to 
feed the output signal into a root-mean— square voltmeter.. The averaging 
characteristics of the root-mean— square meter were checked and, for a 
narrow— band-width signal, were found to give the same mean values as 
those calculated from thermocouple readings. 

-v 

With this method the energy of the u' ^—fluctuation was measured at 
different points across the channel. At = 0.4 the spectrum of the 

turbulent shear u'v' was also obtained. 3 An x— type meter was used, 
both wires making an angle of 45° with the free— stream direction. 

If ej_^ and e ^ are the mean— square voltage fluctuations (corrected 
for time lag) across the two wires it may be shown (reference 12) that 

ej^ — = Constant x u’v* 


If e^_ and &2 are fed into the wave analyzer then the above relation 
holds for all band widths dn; that is. 


eq(n)2 — ep(n)2 = Constant x u'v’ (n) 


^This method appears to have been first used by Dr. Stanley Corrsin 
in a turbulent jet (private communication) . 




NACA TN 2123 


21 


PRELIMINARY MEASUREMENTS IN A 1-INCH CHANNEL 


In the Initial stages of the turbulent-channel- flow investigation 
a two-dimensional channel of 1— inch width was used. Measurements of 
mean fluctuating velocities and of the correlation coefficient k 
were completed. The scales Iy and Lz were also measured at the 
channel center and were found to be about 0.2 to 0.3 centimeter. This 
small-scale turbulence existing in the channel imposed a definite limi- 
tation on the accuracy of the fluctuating measurements. Wire— length 
corrections as high as 30 percent had to be applied to the measurements 
of v' and w' . Furthermore, no microscales could be measured accu- 
rately; thus one of the objectives of the investigation, the calculation 
of the energy dissipation across the channel, could not be obtained. 
Nevertheless, results show good consistency; the three independent 
measurements of the turbulent shear indicate satisfactory agreement. 

It is of interest to present these preliminary measurements, the 
Reynolds number of which was 12,200, and to compare them with those 
obtained in the 5~ inch channel. Figure 5 shows the distribution of all 
the measured quantities in the 1— inch channel and it may be directly 
compared with figure 1 6 which shows the quantities measured in the 
5— inch channel at nearly the same Reynolds number. 

Although the slopes of the mean velocities at the wall are aim 
the same in each case [for the 1— inch channel — = 0.0019; for the 

V pUo 

5— inch channel T ° = 0.0018 ], the velocity ratios u/U D in the 5— inch 

PUo 2 ) 

channel seem to be lower across most of the channel. This indicates that 
a more intense * turbulent-energy production t|^- takes place in the 

5— inch channel. Indeed, the turbulent— velocity fluctuations u’ /U Q , 
v ’ /U 0 , and w' /U G (especially u'/U 0 ) increase at a faster rate toward 
the wall in the wider channel, although their values at the channel 
center check very well in the two cases. The shear distribution t / pU 0 ^ 
is almost the same for the two channels; it follows, therefore, from 
the relation 


u’v* 


Un 


^ U 1 V 1 

U 0 U 0 
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that k must have a higher absolute value in the 1— inch channel since 
both u' /U 0 and v* /bo are smaller there. Indeed the maximum, value 
of k in this case is -0.6, while in the 5— inch channel it is -O.5.- 

As a matter of interest it could be mentioned that if the basis of 
comparison is not Reynolds number but maximum mean speed (for the 1— inch 
channel Uo = 15 m/sec) then figure 5 shows a distribution of u* , v’ , 
and w’ very similar to that of figure 18 , the absolute values being 
somewhat lower in the 5— inch channel. 


RESULTS AND DISCUSSION 
Mean— Velocity Distribution 


A careful study of the two-dimensional nature of the ch ann el flow 
was first made. Mean-velocity measurements were carried out at approxi- 
mately x = — 2 inches, at different heights in the channel: At positions 

6 inches from the bottom, 6 inches from the top, and at the middle. 
Agreement among the three * sets of measurements confirmed the two- 
dimensionality . A further check was made on possible end effects that 
might influence the results. The length of the channel was extended 
by 6 inches and the mean— velocity measurements were repeated at 
x = -8 inches. • No change in the profile was noticed. Figure 6 shows 
the mean— velocity distributions at three Reynolds numbers, R = 6l,600, 

R = 30,800, and R = 12,300. The distributions follow Von K&mi&n's 
logarithmic law very well, except, of course, near the wall and at the 
center of the channel (fig. 7) • The values of U T were obtained from 
the velocity gradients at the wall. 

Measurements were made with special care close to the wall. 
Velocities at a large number of points were recorded within the laminar 
sublayer in order to establish with reasonable accuracy the shape of 
the velocity profile at the wall (fig. 8). The thickness of the laminar 
sublayer (the point where the velocity distribution deviates from the 

logarithmic law) was found to be 5 ~ 30 (fig* 7)* 

Or 

In figure 8 it can be seen that for the case of the lowest Reynolds 
number a few points near the wall indicate very low velocities. Since 

the point at which ^ = 0 has been determined with great accuracy and 

since the equations of motion of the channel flow require a negative 
curvature for the velocity profile, it was concluded that the hot-wire 
indicated too low velocities near the wall for the lowest Reynolds 
number. (The dashed curve for this case in fig. 8 shows the interpolated 
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velocity distribution for 0 < £ < 0.01.) It can be shown that veiy 

high local— velocity fluctuations cause the hotr^wire to read velocities 
lower than the true mean value . Since the mean current of the wire 
varies with the fourth root of the velocity. 



where u is the velocity to which the ho t-wire responds and- u is the 

true mean velocity. Near, the wall ^ < '- < ^ ^ can therefore be 

neglected. Because of the cosine-type directional characteristics of the 
hot-wire the effect of the w* — fluctuation can also be neglected. The 
velocity' measured by the wire can be written, after expanding the above 
expression, as 


u 



+ • 



Since the average of the odd terms is zero the series converges rapidly. 
In the region in question, the velocity fluctuations obtained were very 


high indeed 


( Ay 


> 0.3oj. Their absolute values, however, are probably 


even larger, since the nonlinearity effect of the large fluctuations on 
the hot-wire response still further amplified by the very low mean 
velocities is not taken into account because the wire— length corrections 
are neglected. The mean-velocity correction therefore given by the 
above relation is too small to give the required 'negative curvature of 
the profile. 


Turbulence Levels 

Figures 9 to l4 represent the results of measurements of the three 
components u’ , v’ , and w' of the turbulent-velocity— fluctuation 
distribution in the channel for the three Reynolds numbers. 

The velocity fluctuations u' relative to local speeds increase 
very rapidly near the wall as is shown in figures 9 and 10. Measure- 
ments very close to the wall indicate that u 1 /u reaches a maximum 

within the laminar boundary layer ^yU T /v ~1 t) and it tends toward 

a const an t value at the wall which is independent of the Reynolds number. 
This point will be discussed in detail under Reynolds Number Effect. 
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The absolute values of the distribution of u' show the same 
general shape as that obtained by Reichardt (reference 9) , having the 
characteristic maximum near the wall and thus showing the strong action 
of viscosity even for values of y a* 46. 

Using the x— type hot-wire technique for obtaining the velocity- 
fluctuation components v' and w', no measurements could be obtained 
near the wall. Figures 13 and 14 show that, while in the center of the 
channel the magnitudes of v' and w' are the same, w' increases 
faster toward the wall. This agrees with the ultramicroscope measure- 
ments in a pipe by Fage and Townend (reference 18) . 

No length corrections were necessary to the measurements of u' 
except near the wall; however, no corrections were applied in this 
region since no measurements of \ z could be made. In this region, 
furthermore, the fluctuations are very large and the values given in 
figure 10 must be accepted with reserve. The hot-wire response for 
large velocity fluctuations is not well understood yet and no correction 
was attempted. Length corrections were applied to the measurements 
of v* and w' . » 


Correlation Coefficient and Shear Distribution 

The correlation coefficient is fairly constant across most of the 
channel (fig. 15) as indicated already by Wattendorf and Kuethe (refer- 
ence 6). The maximum values of k obtained decrease slowly with 
increasing Reynolds number and thus show a definite Reynolds number 
dependence. Existing results indicate, however, that parameters other 
than Reynolds numbers also influence the value of k. The following 
table shows magnitudes of the maximum correlation coefficient as obtained 
by different investigators working with channels of various widths and 
with various Reynolds numbers: 


Experiments by 

Channel width 
(cm) 

R 

^max 

Reichardt (reference 9) 

24.6 

8,000 

-0.45 

Laufer (present paper) 

2.5 (1 in.) 

12,200 

-.63 

Laufer (present paper) 

12.7 (5 in.) 

12,300 

-.50 

Wattendorf (reference 7) 

5-0 

' 15,500 

-.52 

Laufer (present paper) 

12.7 (5 in.) 

30,800 

-.45 

Laufer (present paper) 

12.7 (5 in.) 

61,600 

-.40 

Wattendorf obtained his value of 

k from his measurements 

of the quan- 


titles u'/U 0 , v'/u 0 , and u'v'/Uq^. Unfortunately his 



NACA TN 2123 


25 


measurements of v' are incorrect; their magnitude is larger than that 
of u' at all values of y across the channel contrary to results 
obtained by Reichardt and by the present writer. Wattendorf himself 
points out the inaccuracy of the values of v* in his paper. Because 
of the large values of v* his computed correlation coefficients are 
too low. The value ^max = ~ 0*52 listed in the above table was obtained 

7 ^ i 

by using Wattendorf 1 s values u - - ) ■— = 0.001 and — = 0.055 and the 

U 0 ^ U 0 

value = 0.035 = 0-5^ obtained by the author for both the 

1— inch (2.5-cm) and 5~ : inch ( 12. 7— cm) channel at R = 12,200 and 
R = 12,300, respectively. 

The variation of k indicates that for the same Reynolds number 
the correlation coefficient tends to increase with increasing maximum 
mean velocity (i.e., in a channel of decreasing width). 

Figures 1 6 to 18 show the distribution of all the measured quan- 
tities including the shear distribution. For the case of the two 
higher Reynolds numbers, three independent measurements were made for 
the determination of the shear distribution by methods . indicated in 
Equations of Motion for Two-Dimensional Channel Flow. Figure 19 
indicates that consistent results are obtained for the value of T 
whether calculated from the pressure gradient along the channel or 
from the velocity gradient at the wall. The turbulent-shear distri- 
butions obtained from hot-wire measurements and those calculated from 
figure 17 show satisfactory agreement. However, for the flow at the 
highest Reynolds number, the hot-wire measurements gave a 20-percent— 
lower value for the shear coefficient, as may be seen in figure l8.^ 

At present, the writer can give no satisfactory explanation for this 
discrepancy. In the flow at the lowest Reynolds number, no pressure 
measurements were made because the very low pressure gradient 
(approx. 0.0003 mm of alcohol/cm) did not permit accurate measurements. 
However, in this case the wall shear computed from the measured mean 
fluctuating quantities u’ , v' , and k and the mean— velocity 
meaaurements may be compared with that obtained in the 1— inch— channel' . 
flow having nearly the same Reynolds number. The comparison shows 

T n 

satisfactory agreement: For the 1— inch channel — = 0.0019; for 

pUic> 

the 5— inch channel — = 0.0018. 

pTJo 2 


Nlhis 


result was accepted after the absolute values were carefully 
checked; the v’-wire response was compared with u’— fluctuations in an 
isotropic field and the two-dimensional character of the flow was 
checked. 
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Scale and Microscale Measurements 

For a further understanding of the structure of the turbulent 
field, correlations of u’— fluctuations at two different points were 
carried out. Since the field is not isotropic, the scale and microscale 
measurements were repeated both in the y— and z-direction for different 
values of y/d. 

E z — correlations .— Figure 20 shows typical E z — correlation distribu- 
tions at different values of y/d corresponding to E = 30,800. For 
larger values of z, inaccuracies in the measurements did not permit 
the exploration of the negative region of the correlation distribution. 
From the measured distributions of E z at four stations across the 
channel and for different values of E, the values of «.n<i x z 
were calculated (figs. 21 and 22). In these figures 1^ is seen 
to be decreasing uniformly toward the wall, while X z reaches a definite 

y 

maximum at about — = 0.7 and then decreases with decreasing y/d. 

Ey— correlations Distributions of the By— correlation at a given 
■ " ■ ■ 1 - \ 

value of y/d were obtained by fixing the stationary hot-wire at the 

given value of y/d and traversing with the moving hot-wire away from the 

fixed one toward the channel center. Values of Iy and Xy were then 

calculated from these distributions of E y . In the region 1.0 > — > 0.1 

J d 

the gradients of the various quantities are slight; no significant asym- 
metry in the function Ey is therefore expected. Measurements of the 
top part of the ^-correlation distribution (l.O > Ey > 0.8) by Prandtl 
and Beichardt (reference 19) Justify this belief fairly well. 

By comparing figure 20 with figure 23, the strong effect of the 
existing nonisotropy can immediately be seen. Although both Ey. and E z 

correspond to Von Harman* s g— function (reference 10) they exhibit dif- 
ferent behaviors across the channel. The B z -correlation function falls 
more and more rapidly to zero with decreasing y/d. Traversing 

from. ^ = 1.0 to ^ = 0.70 the distribution of Ey is found to behave 
similarly; however, the dashed curve measured at £• = 0.4 (fig. 23) 

indicates considerably increased correlations for larger values of T. 

It is thus seen that around ^ = 0.7 there is a definite decrease in 

energy content of the fluctuations having low frequencies. Further con- 
sideration of this fact is given later in the discussion of the energy 
balance in the channel. 
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n 00 

The values of Iy = j By dY (fig. 21 ) are not so reliable as 

those of Lz since By could not be measured at sufficiently large 
values of T because of the limitation of the. traversing mechanism 
used. 


The distribution of Xy across the channel is similar to that 

of X z (fig. 22). From £ = 1.0 to w 0.7 both Xy- and X z 

d d J 

increase almost proportionally to u* , indicating that the turbulent- 
energy dissipation W ~ is approximately constant in this region. 


It should be noticed, however, that X z is consistently larger than 
throughout the channel cross section. Some measurements of Xy close 
to the wall are also indicated in figure 23 . No length correction was 
found to be necessary for the measured values of Xy and X z . 


A rough est ima te of Lk by the method already described gave a 
value approxi m ately twice that of 1^ at the center of the channel 
(fig. 21). Its value increases to a maximum at 0.5 and then 

decreases rapidly. No values for L*- are given for the lowest Beynolds 
number; in this case the value of a is very close to unity and, since 

Ljr is proportional to — — y , a should be known within an accuracy 

of 1 percent to give consistent results for L*-. Unfortunately measure- 
ments of Ljj- cannot be made with this accuracy, particularly when the 
u'— fluctuations are of rather low frequencies as is the case 
for E = 12,300. It should be mentioned that the accuracy of the deter- 
mination of Ljj- is more limited by the inaccurately measured value 
of a than by the fact that an approximate spectrum function F(n) is 
used in equation ( 9 ). 


The distribution of obtained by the differentiation method 

shows the same behavior as that of Xy and X z (fig. 22). It ha s 

y 

the characteristic maximum at -■ £ 0.7* No measurements were made 

for R = 12,300 "because of the inaccuracy due to the large noise— to— 
signal ratios. 

As a matter of interest the results of some measurements of 

obtained with the zero-counting method are compared with those described 
above. (See fig. 22.) Since the validity of the assumptions involved 
in this method is not clarified yet, these measurements should be 
regarded as preliminary. 
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Spectrum Measurements 

The spectrums of the velocity fluctuations u* 2 have been obtained 
at various values of y/d' across the cha nn el. With Improved instru- 
mentation it was possible to reduce experimental scatter by a consid- 
erable amount. The accuracy of the present measurements is believed to 
be within ±10 percent with the exception of values corresponding to low 
frequencies (n < 100 cps) because of the large-amplitude fluctuations 
and with the exception of values corresponding to high frequencies 
(n > 4000 cps) because of noise and of possible wire— length effect. A 

typical spectrum distribution taken at j = 1.00 is shown in figure 2k. 

(The measured spectrums at various positions of y/d are given in 
table I.) No attempt is made to compare the measured distributions with 
existing theories since the restrictive assumptions of these theories 
(isotropy and flows at high Beynolds numbers) are not satisfied in the 
present experiments. The only frequency range where comparison is 
possibly Justified is the viscous region, that is, the high-frequency 
part of the spectrum where viscosity plays a dominant role. Unfortu- 
nately here the accuracy of the measurements is not good enough to afford 
any definite conclusions. Therefore, the rather close agreement of the 
measured spectrum (for all values of. y/d except in the laminar sub- 
layer) with the n~ law predicted by Heisenberg (reference 2.) should be 
accepted with reserve (fig. 2k). 


For a comparison of the energy spectrums at different points in 
the channel, figure 25 shows the distribution of n 2 F(n) instead of 
that of F(n) as a function of the frequency. This method of presen- 
tation emphasizes energies corresponding to the higher frequencies but 
gives a clearer over-all comparison of the different sets of measure- 
ments. The characteristic maximum of the distribution of imme- 

diately becomes apparent in this figure. From the fact that the 
correlation coefficient is a Fourier transform of the spectrum function 
it follows that 



1 kn 2 r°° 2w \ ^ 
— = -5- / n^F(n) dn 

-4: u £ J 0 


Thus X, x can be computed from the distributions of n2F(n) . The 
calculated values check, within a few percent with direct measurements 
(fig. 22). 
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Reynolds Number Effect 


Figure 7 shows the distribution of mean velocities plotted in 
the form of "friction velocity" u/U T against "frictionr-distan.ee 
parameter" yU T /v. In this form the profiles are independent of the 
Reynolds number U 0 d/v and follow Von Karmen's logarithmic law: 


yUT 

^ = A log - — + B 
Ut V 


The constants A and B are 6.9 and respectively. Comparing 

the present measurements with previous ones it is to be noted that 
because of the very smooth wall surfaces used, the value B is larger 
than that of other investigator's. Making allowances for this effect, 
measurements are in fairly good agreement with those of Doench (refer- 
ence 4) and Reichardt (reference 9)- Nikuradse’s channel results 
(reference 5) differ from all other existing experiments. It is 
pertinent to mention that for the Reynolds number range in question . 

(R < 100,000) the slope of the logarithmic mean— velocity distribution 
(the constant A) both in pipes and channels is larger than that 
established by Nikuradse's pipe measurements (A = 5»75)» 

The Reynolds number has a definite influence on the velocity 
fluctuations also. Along most of the cross section where the influence 
of the viscosity is negligible the fluctuations u', v' , and w' 
decrease slowly with increasing Reynolds number. 


Very close to the wall, according to Prandtl’s hypothesis, all 
velocities of the form Velocity /U t must be a function of the friction- 
distance parameter only. It has already been pointed out that the mean 
velocity obeys this similarity law. Figure 26 shows the distribution 
of u' /U T as a function of yU-r/v for the different flow conditions 
investigated. The following remarks can be made with reference to this 
figure: 


(l) For values, of 


yu T 

— < 100 , 

V 


values of u' /U T corresponding to 


the various cases indicate similar behavior. They reach a maximum at 

yU T 

about % 17. 

V 
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(2) It Is believed that the similarity is actually more complete 
than indicated in figure 2 6 . Because the microscale is not known in 
regions very close to the wall no length corrections could be applied. 
These corrections would of course be appreciably higher for the two 
high-velocity flows (R = 61,600, 2d = 5 in. and R = 12,200, 

2 d .= 1 in.), and would therefore bring the various distributions of u'/U T 
closer together in the region in question. 

( 3 ) The effect of viscosity is more pronounced on the fluctuating 
quantities than on the mean velocity. 

(4) Taylor pointed out in 1932 (supporting his arguments, by Fage 
and Townend's ultramicroscope measurements) that u'/u and w'/u 
approach a finite value at the wall (reference 18). It follows from 
the similarity law that this value should be an absolute constant inde- 
pendent of the Reynolds number. Figure 10 indicates this to be true, 
the constant being (u 1 /u)y_Q ~ 0 .l 8 . 

It is of considerable interest to discuss the variation of the 
scale and microscale with Reynolds number. For flows behind grids where 
the turbulence is isotropic the scale is independent of the mean velocity 
and depends on the mesh size of the grid. Similar behavior was fo und for 
the c hann el flow. Figure 21 shows the distributions of Iy anr) 1 ^. for 

different velocities. These distributions indicate no consistent varia- 
tion with velocity. Furthermore, measurements in a 1— inch channel give 
a value for I^ five times lower and a ratio for Iy somewhat larger 
than the values obtained in the present investigation. 

The variation of X depends, of course, on the velocity and c hanne l 
width. The values of X decrease with increasing velocity; however, the 
variation of X with channel width is less than that of L. 


Fully Developed Character of Turbulence 


The flow in the channel is called fully developed if the varia- 
tions of the mean values of the velocity and the mean squares of the 
velocity fluctuations with x are very small. That the mean velocity 
profile does not vary downstream is evident from the pressure-gradient 

measurements (fig. 19). The gradient in x of u ' 2 was measured on 
the axis of the channel. It was found that u * 2 was indeed decreasing 

with x. The gradient, however, was very small as compared with — 

P dx 


• du ,g 
dx 


R* 0.01 


1 dp 
p dx 
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Hence for all practical purposes can be neglected. No 

measurements have been made concerning du’v’/dx, since the scatter 
in the values would cover any effect. However, there is little doubt 

that chi'v' /8x is of the same order as du’^/dx and that the use of 
equations (3a) and (3b) is therefore justified here. 


Energy Balance in Fluctuating Field 


The energy- equation for a two-dimensional channel has the form, 
given by equation (8) and is valid throughout the cross section of 
the channel with the exception of a small region near the wall. 


The term on the left side of equation 


(8) corresponds to the 


energy produced by the shearing stresses and it can be obtained directly 
from the measurements of t and from the mean velocity profile. The 


second term on the right jil 


/chyWduj'N 

IW j I 


expresses the amount of energy 


that is being dissipated because. of the breaking down of the larger 
eddies to smaller ones. The term may be written explicitly 


W = 



+ + fdw ' , dv* \ 2 + /du* + 8w'\ 2 . 

8y / W 8z / \8z 8x / 


The problem is to express these functions in terms of easily meas- 
urable quantities. In the case of isotropic turbulence Taylor solved the 
problem by introducing the microscale of turbulence X, and obtained for 
the dissipation 

W = 15H (10) 

' 4 

It is attempted here to carry over his analysis to the case of 
channel- flow. The following assumptions have to be made in this 
connection: 

(a) The gradient of the velocity fluctuations is smell in all three 
directions. With the exception of the region near the wall this is 
Justifiable from the measurements. 
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(b) Since only correlations of u' have been measured, it is 
assumed that the correlations of w* as functions of Z (Von Earman’ s 
f— function) have the same curvature at the origin (T = Z = 0) as the 
correlation of. u* has as a function of X; that is, 

(^L\ J^L) ._jl 

\ bZ 2 Jo \ bX 2 Jo 

Furthermore the correlations of v* as functions of X and Z and 
the correlations of w' as functions of X (Von Kar me n 's g-functions) 
have the same curvature at the origin as the correlation of u* has as 
a function Z; that is, 



Finally, 

/;>y’\ /ay\ /ay\ g 

\ St 2 Jo \ bY 2 ' 0 \ bY 2 Jo Xy 2 

(c) The cross products were calculated using similar arguments; 

thus , 

<3v' du‘ _ _ 1 V* u* 

dx dy 2 X z Xy 

dw 1 dv 1 1 w 1 v* 

by dz 2 Xy X z 

c>u» dv’ _ _ 1 u^ 

<3z dx 2 X z X z 
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With these as sumpt ions the derivatives of the fluctuations could 
he expressed in terms of the measured values of and X z : 


W = n 


-i-p + ^5) + 2im + 2 im_iiiL + 2 “l! + a2l£_ii2i + 

y V 2 V V 2 V 2 V. 


v'2 u'2 u , w' w 1 2 

2 - — + 2 — + p - — 

2 2 o ^ ? 

Xj> X z X z X z 


At the middle of the channel W turned out to he 2.86 ergs per cubic 
centimeter per second for E = 30,800. Using the isotropic relation 


W = 




the values 2.88, 4 . 38 , and 2.24 ergs per cubic cent ime ter per second 
were obtained depen ding upon whether X^., Xy, or X z was used. 

(For the value of u'^ the algebraic mean of the squares of the fluc- 
tuations was used.) 

Figure 27 shows the distribution of W in the center region of 
the channel. Taylor obtained a similar distribution of the dissipated 
energy across the channel (reference 20 ) ; however, his numerical magni- 
tudes are too high since he substituted in the isotropic relation, in 
equation (10), the values of Xy, which turn out to be less than those 
of X^ and X z . 

From the known distributions of the energy production and dissi- 
pation the diffusion of energy is easily calculated from equation (8). 
It should be pointed out that because of the approximations involved in 
the calculation of the dissipation and diffusion terms, the conclusions 
derived from them are more or less of a qualitative nature. It is seen 

from figure 27 that at ^ s: 0.7 the diffusion term is zero. It was 

also pointed out earlier that in this region the Ry— correlations show 
a considerable ch an ge in shape indicating a shift in energy from the 
lower to higher frequencies of the velocity fluctuations. These two 
facts suggest that the energy diffusion is associated mainly with the 
low frequencies of the fluctuations. 
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The equation expressing the balance of the three forms of energy 
furnishes the following picture of the turbulent flow field in the 
channel where viscous dissipation is still negligible: Two planes 

passing through points where the diffusion of energy vanishes divide 
the channel flow into three parts. From these planes energy is being 
transported toward the channel center and the walls. The middle region, 
the width of which is of the order of I^, receives most of its energy 
by diffusive action and this energy is dissipated here at a constant 
rate. In the two outside regions all three energy terms increase 
rapidly, the production term being the dominant one, and their inter- 
action is more involved. 

This picture of the flow field is only of a descriptive nature. 

The purpose of further investigations should be to obtain information 
on the development and mechanism of such an energy balance. 

Locally Isotropic Character of Turbulent Channel Flow 

The concept of locally isotropic turbulence obtained by Kolmogoroff 
requires the smaller eddies in turbulent flow to approach isotropy. 
Small er eddies are the oneB with length dimension l small compared 
with the scale of turbulence L. The smallest characteristic length in 
turbulent motion is Kolmogoroff 1 s q defined as 



Where e is the total dissipated energy. Clearly, to approach locally 
isotropic conditions it is necessary that 

L » q 

It -follows from this hypothesis . that beyond stuff iciently high 
frequencies (of the order of u* /q, say) no correlations exist between 
the components of the velocity fluctuations. Figure 28 shows the meas- 
ured spectrum of u' 2 as compared with the spectrum of u’v’ 

at 21 = o.JU-. It is seen that the shear spectrum tends to zero at a 
d 

frequency of about 1500 cycles per second while F—yg (n) still has a 

detectable value at 5000 cycles per second. This result verifies 
Kolmogoroff 1 s assumption. It should be mentioned that the absolute 

values of F (n) are not so accurate as those of F— 75 (n) since 

u'v 1 u ,£ = 

they represent the small difference between two large values of hot- 
wire signal. 
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It is to be noted that even though local isotropy was shown to 
exist i n the channel flow, the values of the various vorticity 


terms 



differ appreciably, indicating that the relative magni- 


tudes of these terms do not constitute a sensitive test for the existence 
of local isotropy. This is evident from the comparison of the distribu— 



(fig. 28 ). It is apparent from the 


figures that a large . part of the contribution to 



comes 


from a frequency well below 1500 cycles per second. Local isotropy, 
according to the shea? spectrum, exists only above a frequency of 
1500 cycles per second. 


CONCLUDING REMARKS 


The measurements presented here confirm the general conceptions 
concerning the mean velocity profile in a turbulent channel. The extent 
of the laminar sublayer, the velocity profile in the sublayer, and the 
over-all velocity distribution as measured here are in good agreement 
with general theoretical expectations. 


Measurements of the turbulent field show that the hot-wire tech- 
nique is well enough developed to give consistent results for the 
intensities and correlation functions. 


Detailed measurements of the velocity fluctuations in the direction 
of the flow u' were carried out well within the laminar sublayer. It 
was found that the similarity law for u'/Ut, where Ut is friction 
velocity, holds fairly well^in the vicinity of the wall and as a con- 
sequence the magnitude of u'/u, where u is local mean velocity, 
approaches an absolute constant at the wall (approx. 0.l8), the value 
being independent of the Reynolds number. The magnitudes of the 
velocity fluctuations normal to the flow in the lateral and vertical 
directions v' and w' are nearly the same in the middle region of 
the channel, w' increasing more rapidly toward the wall. 


The measured microscales of turbulence X x , Xy, and X z con- 
sistently show a maximum at % 0.7, where y is lateral distance 

and d is half width of the channel; they increase proportionally 
with u’ , indicating a constant rate of ene rgy dissipation W in the 

u* 2 

center portion of the channel since V ~ 
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The scales of turbulence Iy and in the center region are 

independent of Eeynolds number and depend only on the channel width. 

The microscales, however, show a dependency on the Eeynolds number. 

From the calculated magnitudes of the energy produced by the 
shearing stresses and of the dissipated energy a descriptive picture of 
the energy diffusion in the center region, of the channel is obtained. 

The spectrum measurements of the u'^— fluctuations tend to indicate 
that (n) behaves as n - ? over a large band in the high-frequency 



•b and width dn and n is frequency). From the comparison of the 

spectrums of u’ ^ and of the turbulent shear the existence of local 
isotropy in the channel flow is verified. 


Guggenheim Aeronautical Laboratory 

California Institute of Technology 

Pasadena,. Calif., September 1, 19^9 
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TABLE I.- MEASURED VALUES OF THE SPECTRUM OF u ' 2 
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( sec) 
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250 
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Figure 2.- 


NACA. 

Exit of channel. 


Page intentionally left blank 


Page intentionally left blank 




Figure 3.- 


..NACA. 

Traversing mechanism. 



Figure 4.- Typical oscillogram of u ' -fluctuation at 
y/d = 1.0. Horizontal line corresponds to u'(t) 
approximately 1/20 second. 
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Figure 5.- Distribution of measured quantities across two-dimensional channel of 1-inch width. R = 12,200. 












kQ 



















measured 


'> M o 
0 & 



f 




[ 


n 

i*lo° 




Y 

>|p° 



Igure 16.- Distribution of mean and fluctuating quantities. 





calculated from velocity distribution 



Figure 17.- Distribution of mean and fluctuating q uant ities. 
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1.55 centimeters, = 0.81 centimeter. 
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Figure 26.- Logarithmic representation of velocity fluctuations 




Figure 27.- Turbulent-energy balance in center region of channel. 






Figure 28.- Comparison of spectrums of u 1 ^ and uV at y/d 



